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Abstract 

The superconductivity of quasi-one-dimensional electrons in the magnetic 
field is studied. The system is described as the one-dimensional electrons 
with no frustration due to the magnetic field. The interaction is assumed 
to be attractive between electrons in the nearest chains, which corresponds 
to the lines of nodes of the energy gap in the absence of the magnetic field. 
The effective interaction depends on the magnetic field and the transverse 
momentum. As the magnetic field becomes strong, the transition temperature 
of the spin-triplet superconductivity oscillates, while that of the spin-singlet 
increases monotonically. 
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running head: Superconductivity in strong magnetic field 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The reentrance of the superconductivity in quasi-one- dimensional electrons in strong 
magnetic field attracts theoretical interest [ffl-B]. In the quasi-one-dimensional systems the 
Fermi surface is open and Landau level quantization does not exist. Lebed' has shown that 
the superconducting transition temperature reaches to that in the absence of the magnetic 
field as the magnetic field is increased. Dupuis et al. have shown that there exists the 
cascade of the superconducting phases separated by first order transitions. The reentrance 
of the superconductivity is understood as follows. When the strong magnetic field is applied 
along the b axis, the semi-classical orbits of electrons are localized in the a-b plane, where 
a is the most conducting axis, i.e. the motion along the c axis is bounded. Since the 
orbital frustration due to the magnetic field comes from the motion of electrons in the plane 
perpendicular to the magnetic field, electrons can make Cooper pairs without affected by 
the orbital frustration and the transition temperature increases. 

The reentrance of the superconductivity in two-dimensional electrons in strong magnetic 
field also attracts theoretical interest 0,0. This phenomena is similar to that in quasi- 
one-dimensional systems, but the origin for the reentrance of the superconductivity in two- 
dimension is different from that in quasi-one-dimension. In two-dimensional case Cooper 
pairs in strong magnetic field are formed between electrons at the lowest Landau level, 
whereas superconductivity is stabilized in magnetic field without Landau level quantization 
in quasi-one-dimensional system. 

Quasi-one- dimensional electrons are realized in organic conductors, (TMTSF^X, where 
anion X is PF 6 , CIO4 etc. The system can be described as the anisotropic tight-binding 
model with the hopping matrix elements, t a ~ 0.25eV, tb/t a ~ 0.1 and t c /t a ~ 0.01. Due to 
the smallness of t c , the magnetic field necessary to the reentrance of the superconductivity 
is estimated to be about 20T Jl]||, which can be realized in experiments. 

The superconductivity in the quasi-one-dimensional organics, (TMTSF) 2 X, is thought 
to be not the conventional BCS type even in the absence of magnetic field. When the nest- 
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ing condition of the Fermi surface is satisfied sufficiently, which is thought to be the case 
at low pressure, the quasi-one- dimensional organic conductors are not superconductors but 
insulators due to spin density wave (SDW). Since SDW is caused by the on-site Coulomb 
repulsion, the on-site repulsion is expected to be large in this system. As a result the con- 
ventional s-wave spin-singlet superconductivity is not likely to be stabilized. Indeed, the 
NMR relaxation rate, T^ 1 of the (TMTSF) 2 C104 shows little enhancement just below the 
transition temperature and it decreases as T 3 as temperature becomes low, which strongly 
suggests the line nodes of the energy gap in the superconducting state |§J| . The line nodes 
in the superconducting states are caused by the attractive interaction between electrons 
which are not on the same site. When the magnetic field is applied along the y direction, 
the k y dependence of the interaction is not affected by the magnetic field. Thus if the attrac- 
tive interaction depends only on k y , the magnetic field dependence of the superconducting 
transition temperature is the same as that in the case of the on-site attraction. However, 
the k z dependence of the interaction makes the different field dependence. In the quasi-one 
dimensional case the k x dependence of the interaction is not very important because Fermi 
surface consists of two planes, k x « ±kp. Therefore, in this paper we assume the attractive 
interaction between electrons in the nearest chains in z direction, which results in the line 
node of the energy gap in the superconducting state. 

In the previous paper one of the authors has shown that if the magnetic field is ap- 
plied, the quasi-one-dimensional electrons in the perpendicular plane is described as the 
one-dimensional system and the effective interaction depends on the magnetic field and the 
transverse component of the wave vector [IIJ. When the magnetic field is applied along 



the c axis, the perfect nesting is recovered as long as we neglect the small imperfectness 
of the order of t\jt a . In that paper the on-site repulsive interaction is assumed and the 
field-induced spin density wave (FISDW) is studied. The quantum Hall effect in FISDW is 
understood as a result of the /c^-dependent phase of the effective interaction. 

In this paper we apply the same procedure to the case of the coexistence of the on- 
site repulsion and the nearest-site attraction. When the magnetic field is applied in the y 



direction, the imperfectness of the nesting of the Fermi surface is not changed and FISDW 
is not stabilized. In that case the superconductivity is caused by the attractive part of the 
interaction. 

II. QUASI-ONE-DIMENSIONAL ELECTRONS IN MAGNETIC FIELD 

We study the quasi-one- dimensional electrons in the magnetic field. For simplicity, we 
take a, b, and c axes to be perpendicular each other and to be along the x, y and z directions, 
respectively. The generalization to the non-orthogonal lattice can be performed as in the 
FISDW case [|TTJ. The Hamiltonian is written as 



n 


= T~Co- 




Ho 


= ~t a 






-h 












-t c 





MU= E U ij C i,aCi,*Cj,a>Cjy (!) 
<i,j>,CT,cr' 

where c| CT and Cj )(7 are creation and annihilation operators of electrons, t a , tj, and t c are 
the hopping matrix elements along the a, b and c axes, respectively, Uij is the interaction 
between electrons at i and j sites and 

% = ^ [' Adl (2) 

In the above A is the vector potential and 0o = hco/e is the flux quantum, where Cq is 
the light velocity. We consider the anisotropic system with the hopping matrix elements 
t a 3> t b ^> t c , which is the case in (TMTSF) 2 X. The nesting of the Fermi surface is not 
perfect due to the tb term. 

In eq.(0) we have neglected the Zeeman term for simplicity. The Zeeman term does not 
play any important role for the equal-spin-pairing case of the spin triplet. If the Zeeman 
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energy is taken into account, the transition temperature of the spin singlet is reduced but the 
superconductivity is not completely destroyed, since half of the density of states is available 
to make Cooper pairs for the Larkin-Ovchinikov-Frude-Ferrell state |3],^[]. Therefore, the 
effect of the Zeeman energy for the spin singlet can be taken into account by putting the 
density of states to be half in the strong magnetic field. 

In this paper the magnetic field H is applied in the y direction and the vector potential 
A is taken as A = (0, 0, —Hx). Then the non-interacting Hamiltonian is written as 



'He 



2 dkv 
( 



kF -% G 



dky 
2n/b 



E 



i(k) 



dk z 



M_i T 
T* M T 
T* Mi 



T 



Cff(k), 



(3) 



where 



M n = — 2t a {cos[a(k x + nG)] — cos(ak F )} — 2t b cos(bk 
T = —t c exp(ick z ), 

4 = (..., 4(k-G),4(k), C t(k+G)...), 



(4) 
(5) 
(6) 



and G = (G, 0, 0) = {ecH/{hco), 0, 0). 

The creation operators of electrons can be written in terms of the creation operators of 
the eigenstates of the non-interacting Hamiltonian in the presence of the magnetic field. We 
consider the case that the system has the open Fermi surface and the magnetic field is not 
very strong in the sense that the magnetic flux per plaquette is much smaller than the flux 
quantum. Then we get G kp and we can treat the electrons with k x {^ hp) and —k x 
independently in TC as linear combinations of the eigenstates ^l(n, k) and \l/|(n, k) as |TD 
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: f (k + mG) = e imck * Prim, n)^ r {n, k) 



zU-k + mG) 



-imckz 



^2<Pi(m,n)^l(n, -k), 



where m and n are integers. 

In this paper we take the interaction as 



U if i?i = rj 

U\ if r i = Tj ± cz 

otherwise 



The Fourier-transform of the interaction is obtained as 

U(k, k') = U + 2Ui cos[c(k z - k' z )}. 
The interaction Hamiltonian is written as 

nv=Y. EE 

k,k',q m,m' ,N a,a' 

[U(k z , k' z - q z )cl(k + mG)4,(-k + q - mG + NG) 
xc CT /(k' + m'G)c CT (-k / + q - m'G + NG) 

+U(k z , *4)4(k + mG)4,(-k + q - mG + NG) 
xc a t{-\s! + q - m'G + NG)c a {k' + m'G) 



(7) 



(8) 



(9) 



(10) 



where m, m! and iV are integers and kp — G/2 < k x < kp + G/2. The first and the second 
terms are the generalization of the g\ and gi terms in the g-ology of the one-dimensional 
system [|i"2l , |l3| to the quasi-one-dimensional system. 

In the previous paper |TI| particle- hole channel is taken into account, since these terms 
make the instability of the Fermi surface into the spin- density- wave state if the magnetic 
field is applied along the z direction. Here we consider the particle-particle channel, which 
corresponds to the instability toward the superconductivity in the case of the attractive 
interaction. 
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With the eigenstates ^(n, k) and ^i(n, k) the interaction Hamiltonian is written as 



Hv= E E E ^ 

,0-2,0-3,0-4 (ni , h*2> "^3; k z , k z , q z 

k,k',q«i,n2,«3,»M 01,02,03,04 

xtt Ja3 (n 3 , -k' + q)^ rCT4 (n 4 ,k') 



(11) 



where 



E7 



01,02,03 



,a A {ni,n 2 ,n 3 ,n4,k z ,k' z ,q z ) 



i \2mc(k z - |) - 2m'c(A;; - |) - Nc(k z - k' z )} 



= E ex p 

m,m' ,N 

X(p r (m, ni)(pi(—m + N, n 2 ) 
Xip*(—m' + N, n 3 )ipl(m', n 4 ) 



-U{k z ,-k' z + q z )5 

01,03^02,04 4" U{k Zl k z )5 (Tl ^ (T4 5 a2 .,T 



(12) 

The coefficients (p r (m,n) and ipi(m,n) can be calculated by diagonalizing the matrix in 
eq.(|3]) numerically as we have done in the previous paper fllQf . In this paper, instead of 
diagonalizing the matrix numerically, we use the approximation that the dispersion in k x is 
taken to be linear, i.e., 



M n ~ vp{k x + nG — k F ) — 2t b cos(bk y ) , 



(13) 



where Vp = 2t a a sin akp is the Fermi velocity. With this approximation the eigenvalues are 
obtained as 



£(n, k x , k y ) = Vp(k x + nG — kp) — 2tbCos(bk y 



and the eigenstates are given with the coefficients 



(14) 



ip r (m,n) = J-n+m (z) 



(15) 



and 



ipi(m,n) = J n - m (z) 



(16) 



where J n {z) is the Bessel function and z = 2t c /(vpG). We can perform the m and ml 
summation in eq.(|T2D by using the identity, 



e imQ J m _ N (z)J m (z)=e l - {Q+7r) J N (2zsm 



(17) 



m=— oo 



and get 



,0-2,0-3,(74 (^1) n 2t n 3i n 4i k z , k z: Qz) 



J^exp i i (m - n 2 )c(k z - -^) 



+ (n 3 -n 4 )c(A£ - y) 

— (m + n 2 + n 3 + n 4 - 2N) 



71 



xJ N _ m _ m (2z sin c(k z - y)) 
x J N ^ n3 _ ni (-2 z sin c(k' z - —)) 



-U (k z , k z + q z )5 ai :< T 3 fia 2 ,04 4- £^ (^2) ^2)^0-1,0-4^ 



02,03 



(18) 

In the strong magnetic field the instability of the Fermi surface toward forming Cooper 
pairs occurs due to terms with q x — 0, n\ — —n 2 (= n) and n 4 = — n 3 (= n'). Taking only 
these terms corresponds to the quantum limit approximation ||. In this paper we take this 
approximation, which is justified when T <C vpG. 
The order parameter is defined as [114 



^02,01 (p>) k z , q z ) — ,02,03,04(^5 n, n ,n , k z , k z , q z ) 

n',k' 01,02,03,04 

x < ^ la3 (-n', -k' + q)* r<74 (n', k') > . 



(19) 



The mean-field Hamiltonian is written as 



Hmf= J2 [t(n,k x ,k y ){¥ ra (n,k)ty r(T (n,k) 

n,k,o 

+*k(-n,-k)¥ I(T (-n,-k)}] 

- [^a 2 ^,k z ,q z )^t ai {n,k)^l 2 {~n,-k + q) + H.c. 



(20) 



n,k,q,oi,02 



Note that £(n, k x , k y ) does not depend on k z , which means the one-dimensional motion in 
the x-z plane in the semi-classical picture. 



III. TRANSITION TEMPERATURE 



The transition temperature is given by the linearized gap equation, 

f kF+ % d K ft d K f° d K 



tanh 



t(.n',k' x ,k' y ) 



U <n,<T2m,°A n i ~ n i ~ n 'i n 'i k ^ Ki Qz) ocf , J T u s A g3ff4 (w / , k' z , q z ). (21) 
This equation is simplified by defining A m crio - 2 by 

oo 

A a2(71 (n,k z ,q z ) = e 2 ^-^ £ e imc ^"f ) A m , (T2(T1 , (22) 

m=— oo 

and using 

v z-^+f ^ /-f rffc; tagh jgggj = iv(o) 2t^ 

G L f 2-K/b2^n',k' x ,k' y ) 2 11 ttT ' 1 j 

where iV(0) is the density of states for one spin at the Fermi energy, uj c is the cut-off energy, 
and 7 is the exponential of the Euler constant. The factor 1/2 in the right hand side of 
eq.fl2"B|) comes from the fact that only half of the Fermi surface (k x « +kp) is taken in the 
summation. 

A. spin singlet 

We can treat the spin singlet case and the spin triplet case separately. In this section we 
consider the spin singlet case. For the spin singlet case the order parameter is defined by 

A« = \{^ m ,n - A m ,it) (24) 
and the transition temperature is determined by the equation, 

A^ = AT(0)ln^ £ <U« (25) 



m — — oo 
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where 



jyj-(s) __ j c dk z r c dk' z ^-i c (rnk z -m'k' z ) 



m,m' 



2n/c J-?- 2n/c 

C' C' 

x Jo(2z{sm ck z — sin ck' z }) 

x [U + 2U X cos ck z cos c/^] (26) 



If the maximum eigenvalue of the matrix M^f' ,, which we define g( a \ is positive, the tran- 



sition temperature is given by 



r« = 1A o MceM __L_y (27) 



Using the identity for the Bessel function 

r*r d6 



J N (z) = — exp[i(N8-zsm6)}, (28) 
Jo 2n 



we get 



^mm' — —UoL miTn t — UiCm^m', (29) 



where 



r 2n d6 

L m , m > = / — J m {2zsm8)J m ,(2zsm8) (30) 
jo In 



and 



Cm,ra' ^ (-^m+l,m'+l -^m— l,m' — 1 -^m+l,m' — 1 -^m— l,m'+l 

). (31) 

First we consider the case Uq < and U\ — ( i.e. on-site attraction), which has been 
studied by Lebed' WM and Dupuis et al. Noting that 



L m ,m' — -^TO,JV-^Ar,m' (32) 



N=-oo 



where 



T^ = / — e tJW J m ,(2^sin0') 



j\^ff m '^J N {2zs\nck z ) (33) 



2tt/c 



10 



we can show that all eigenvalues of the matrix L m ^ m i are positive or zero. Since L m ^ m i = if 
one of m and m' is odd and the other is even, the matrix (M m m /) is divided into two parts, 
even and odd. When the magnetic field is changed, the maximum eigenvalue is obtained by 
the even part or the odd part, resulting in the cascade transition into the superconducting 
state, as shown by Dupuis et al. In Fig.l we plot the effective coupling constant 

9 /\Uq\ as a function of h = H/H , where H = 2t c hco/ (vpebc). When we take the 
parameters as t a « 3000K, t c « 20K, a « 7A, c 14 A and afcp = 7r/4, we get H « 10T. 
Next we consider the case Uo > and U± < 0. Since 

L m ,m' — ( — l) m -^-m,m' = ( — 1)™ -^m,-m' (34) 

and 

Cm,m' = ( 1) C— ro,m' = ( 1) Cm,-m') (35) 

we get LC = CL = 0, i.e. the eigenstates of M are the simultaneous eigenstates of L and 
C and at least one of the eigenvalues of L and C for each eigenstate is zero. All eigenvalues 
of C m ^ m i are positive or zero, which can be shown by noting 

Cm,m' = ^2 y(Tm+l,N + ^m-l,iv)(^V 1 m+l + ?V,m-l)- (36) 
N=-oo / 

Thus (?^ s ^ is obtained by calculating the maximum eigenvalue of C mim i in the case of t/o > 
and i7i < 0. In this case the transition temperature does not depend on Uq. In Fig. 2 /\Ui\ 
is plotted as a function of /i. In contrast to Fig.l the effective coupling constant does not 
oscillate as the magnetic field is increased. Only in the small region of h ~ 0.2 the largest 
eigenvalues for the even part and the odd part of the matrix C mim i intersect. 

B. spin triplet 

In this section we consider the spin triplet pairing case. The linearized gap equation for 
the spin triplet pairing is given as 
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A®=N(0) \n 1A3uJc Y M (i) ,A«, 



m'=— oo 



(37) 



where A$ is either A, 



■rn 



,TT> ( A ™,U + A m,|T)/ 2 or A m,U and 





(38) 



If the maximum eigenvalue of the matrix , is positive, we define it as and we get 
the transition temperature as 



In Fig. 3 we plot g®/\Ui\ vs. h. The effective coupling constant is given by the eigenvalue 
of the even part or odd part of the matrix , depending on h, as in the case of the spin 
singlet for the on-site attractive interaction. 



In this paper we have studied the transition to superconductivity in the quasi-one- 
dimensional systems in the magnetic field applied along the y direction. The magnetic 
field per plaquette is much smaller than the flux quantum but strong enough so that we can 
apply the quantum limit approximation, which is realized in quasi-one-dimensional organic 
superconductors in the magnetic field of the order of 10T. With this condition the system is 
treated as one-dimension with the effective interaction depending on the magnetic field. We 
studied the case that the interaction between electrons is repulsive for the on-site electrons 
and attractive for electrons at the nearest sites in the z direction. We obtained that the 
effective coupling constant for the spin triplet oscillates as the magnetic field is increased, 
while that for the spin singlet increases monotonically. 




(39) 



IV. CONCLUSION 
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FIGURES 

FIG. 1. The effective coupling constant for the spin singlet in the case of on-site interaction 
(Uq < 0) as a function of the magnetic field h = H/Hq. The effective interaction is given by the 
eigenvalues for the even part of the matrix (solid lines) and the odd part (broken lines). 

FIG. 2. The effective coupling constant for the spin singlet in the case of attraction interaction 
between nearest sites {U\ < 0) as a function of the magnetic field. The solid lines show the largest 
and the second largest eigenvalues for the even part of the matrix and the broken lines show the 
eigenvalues for the odd part. 

FIG. 3. The effective coupling constant for the spin triplet in the case of attraction interaction 
between nearest sites {U\ < 0) as a function of the magnetic field. The solid lines show the largest 
and the second largest eigenvalues for the even part of the matrix and the broken lines show the 
eigenvalues for the odd part. 
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